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Resonant-bar detectors of gravitational wave as possible probe of the noncommutative 

structure of space 
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We report the plausibility of using quantum mechanical transitions, induced by the combined effect 
of Gravitational wave (GW) and noncommutative (NC) structure of space, among the states of a 
2-dimensional harmonic oscillator, to probe the spatial NC geometry. The phonon modes excited by 
the passing GW within the resonant bar-detector are formally identical to forced harmonic oscillator 
and they represent a length variation of roughly the same order of magnitude as the characteristic 
length-scale of spatial noncommutativity estimated from the phenomenological upper bound of the 
NC parameter. This motivates our present work. We employ various GW wave-forms that are 
typically expected from possible astronomical sources. We find that the transition probablities are 
quite sensitive to the nature of polarization of the GW. We also elaborate on the particular type of 
sources of GW, radiation from which can induce transitions that can be used as effective probe of 
the spatial noncommutative structure. 
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Gravitational waves (GW(s)) are small ripples in the 
fabric of spacetime. The present day GW detectors pri¬ 
marily consists of ground-based (LIGO, VIRGO, GEO, 
TAMA etc.) and space-based (LISA) interferometers^. 
However, the search of GWs began with resonant- 
mass detectors, pioneered by Weber in the 60’s2r— . In 
the decades that followed, the sensitivity of resonant- 
mass detectors have improved considerably, though it 
is clear that it could only detect relatively strong sig¬ 
nals from within our Galaxy or the immediate galactic 
neighbourhood^. Nevertheless, the study of resonant-bar 
detectors is fundamental because it focuses on how GW 
interacts with elastic matter causing vibrations with am¬ 
plitudes many order smaller than the size of a nucleus. In 
a bar detector it is possible to detect these tiny vibrations 
corresponding to just a few tens of phonons^, and varia¬ 
tions A L of the bar-length L ~ lrn, with Ah ^ 10~ 19 . 

Interestingly, it has long been suggested in vari¬ 
ous Gedenken experiments that a sharp localization of 
events in space would induce an uncertainty in spatial 
coordinates^ at the quantum level. This uncertainty 
can be realized by imposing the NC Heisenberg algebra 
on the operators representing phase-space variables 

\X {. pj ] — ihSij , [4',;, — iOij — iQC-ij , [p?': Vj\ — 0 ■( 1) 

where 9ij is the constant antisymmetric tensor, which is 
writen in terms of the constant NC parameter 6 and the 
totally antisymmetric tensor tij. Such granularity in spa¬ 
tial structure have been motivated by string theoretic^ 
and quantum gravityiS results also. 

A wide range of theories, dubbed the NC theories, 
have been constructed in this framework. This includs 
NC quantum mechanics (NCQM)ii~— , NC quantum field 
theory (NCQFT)2T— , gauge theories^ and gravity^—. 
Certain possible phenomenological consequences^— 
have also been predicted. Naturally, a part of the en¬ 
deavour is spent in finding the order of the NC parame¬ 
ter and exploring its connection with observations ^ 39 ' 42 . 


The upperbound on the coordinate commutator |0| found 
in2£ is < (lOTeV) -2 which corresponds to 4 x 10 _4o m 2 
for ?i=c=l.— However, recent studies in NCQM revealed 
that the NC parameter associated with different particles 
may not be the same^ii. and this bound could be as high 
as 6 < (4GeV) -2 — (30MeV) -2 —. These upperbounds 
correspond to the length scale ~ 10 _2O m — 10 -17 m which 
overlaps the length scale where the first GW detection is 
expected. Thus a good possibility of detecting the NC 
structure of spacetime would be in the GW detection ex¬ 
periments as it may as well pick up the NC signature of 
space-time. For this purpose we need NCQM of GW- 
matter interaction that can anticiapte the NC effects in 
GW detection events. 

With this motivation, we have studied the interaction 
of GW(s) with simple matter systems in a NCQM frame¬ 
work in^2r^!i. Our particular interest is in the NCQM 
of harmonic oscillator (HO) interacting with GW be¬ 
cause the response of a bar-detector to GW can be cast 
as phonon mode excitations formally identical to forced 
HQ!*. Ii>^ using system parameters typical to the reso¬ 
nant bar detectors and existing upper-bound of the spa¬ 
tial NC parameter we have shown that our NCQM anal¬ 
ysis is likely to be relevant for identifying NC effects in 
the detector read-outs. Thus, NCQM of the HO interact¬ 
ing with GW is of fundamental importance. Therefore, 
we investigate the transition probabilities between the 
ground state and the excited states of this system in the 
present paper by treating the combined effect of GW and 
spatial noncommutativity as time-dependent perturba¬ 
tions. We employ a number of different GW wave-forms 
that are typically expected from runaway astronomical 
events. 

To proceed we first obtain the classical Hamilto¬ 
nian appropriate for the GW-HO interaction system. 
This can be simply done by noting that in the proper 
detector frame the geodesic deviation equation for a 
2—dimensional harmonic oscillator of mass m and fre- 
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quency w subject to linearized GW becomes^ 

mx-i = —mR-’o,kox k — mm 2 x J (2) 

where dot denotes derivative with respect to the coordi¬ 
nate time of the proper detector frame®*, x J is the proper 
distance of the pendulum from the origin and R^o,ko are 
the relevant components of the curvature tensor in terms 
of the metric perturbation h^ defined by££ 

Q\iv — -\- h^ii, |1 << 1 (3) 

on the flat Minkowski background 

The gauge-choice (h 0/J = 0, = 0 ,h£ = 0) that re¬ 

moves the unphysical degrees of freedom (DOF) and ren¬ 
ders the GW transverse and traceless, has been made, so 
that only non-trivial components of the curvature ten¬ 
sor Ri o,fco = —hjk /2 appear in eq.0. The two physical 
DOF brought out by this choice are referred as the x 
and + polarizations of GW and considering 2 to be the 
propagation direction, the surviving components of the 
2x2 matrix hjk in the transverse plane, h\\ = —/122 and 
h \2 = h 21 represent these states respectively. 

Note that, eq. 0 can be used as long as the spacial 
velocities involved are non-relativistic and |ad| is much 
smaller than the reduced wavelength A- of GW. These 
conditions are collectively referred as the small-velocity 
and long wavelength limit and met by resonant bar- 
detectors and earth bound interferometric detectors^ 
with the origin of the coordinate system centered at the 
detector. This also ensures that in a plane-wave expan¬ 
sion of GW, h jk = f (A jk e ikx + A* jk e~ ikx )d 3 k/ (2tt) 3 , the 

spatial part e lkx ~ 1 all over the detector site. Thus 
our only concern is the time-dependent part of the GW. 
If the polarization information contained in Ajk is ex¬ 
pressed in terms of the Pauli spin matrices, hjk takes the 
most convenient form 


Once we have the classical Hamiltonian we can have the 
NCQM description of the system simply by elevating 
the phase-space variables (ad ,pj) to operators ( x\pj ) 
and imposing the NC Heisenberg algebra 0. Note that 
this algebra can be mappe d 35 ' 42 to the standard (9 = 0) 
Heisenberg algebra spanned by the operators X t and Pj 
of the ordinary QM through 

Xi — Xj 2 ^ ^Cij ' Pi — Pi ' ( 7 ) 

Writing the NCQM version of eq.0 and employing the 
map 0 we obtain^ 

p .2 1 , mm 2 

H = + - mzu 2 Xj 2 + T J ok XjP k -^j—9ej rn X J P m 
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_ 2 J^^m^mPk^Qk = ^0 + Pint • (8) 

This Hamiltonian gives the commutative equivalent de¬ 
scription of the noncommutative system 0 in terms of 
the operators Xj and Pj. Since they admit the standard 
Heisenberg algebra, the rules of ordinary QM applies to 
0. Also note that it has been demonstrated in vari¬ 
ous formulations of NC general relativity 24- -- that any 
NC correction in the gravity sector is second order in the 
NC parameter. Therefore, in a first order theory in NC 
space, the GW remains unaltered by NC effects. The 
first two terms in eq.0 represent the unperturbed HO 
Hamiltonian Hq. Rest of the terms are small compared 
to Hq and can be treated as perturbations H ln t- A term 
quadratic in T has been neglected in eq. @ since we deal 
with linearized gravity. 

Defining raising and lowering operators in terms of the 
oscillator frequency vj 



hjk (t) = 2/ (e x a] k + £+a%) (4) 

where 2/ is the amplitude of the GW and (e x ,e + ) are 
the two possible polarization states of the GW satisfying 
the condition e\ +£?_ = 1 for all t. For linearly polarized 
GW the frequency Q is contained in the time-dependent 
amplitude 2 f(t) whereas for circularly polarized GW the 
time-dependent polarization states (e x (t) , £+ (t)) con¬ 
tains the frequency Cl. 

The Lagrangian for the system ©, can be written, 
upto a total derivative term as 

C = —mx 2 — mT^okXjX k — —mm 2 (xj) 2 ( 5 ) 

dj^j 

where R J o,ko =-^ = ~hjk/ 2. Computing the canon¬ 

ical momentum pj = mij — mF J 0k x k corresponding to Xj 
we write the Hamitonian 

H = ^ (Pj + mr o k xk ) + ( x jf ■ (6) 


we write the time-dependent interaction part of the 
Hamiltonian 0 as 


HL(t) = - 


ih ; 

1 
mmO 


hjk(t ) (ajttk - a}4) 


e jmhjk(t) (a m a k - a m a\ + C.Cj 


( 10 ) 


where C.C means complex conjugate. We can now apply 
the time-dependent perturbation theory to compute the 
probability of transition between the ground state |0, 0) 
and the excited states of the 2 -d harmonic oscillator. To 
the lowest order approximation the probability amplitude 
of transition from an initial state |z) to a final state |/), 
(*//), due to a perturbation V(t) is given by^£ 


Ci->f(t -too) = -- 


pt —^-|-oo 


dt' 


Fj k {t')e^ E t- E ^' 


x(^ f \Q jk \^) (11) 
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where V(t) = Fjk(t)Qjk■ Using the above result, we find 
that the probability of transition survives only between 
the ground state 10, 0) and the second excited state and it 
reads 


Co—s-2 T" 

n 



Fjk (t)e^ 2 - E “)‘((2,0|Q, fc |0,0) 

+(l,l|Q jfe |0,0) + (0,2|4 fc |0,0))' 

( 12 ) 


where Fjk it) = hjk{t) contains the explicit time depen¬ 
dence of H' nt and 


Qjk — 


ih 


+- 


mmO 


(ajdk - 

€jm (^m^k “f~ C.C ^ (13) 


contains the raising and lowering operators appearing in 
eq. m- Expanding out Q for i,j = 1,2, we obtain the 
transition amplitude between the ground state |0,0) and 
the second excited state to be 


/• + OO 


Co—s-2 — ~T 


dte 


”“(§W*) + ^ 


hu{t) 


(14) 


The above equation is the main working formula in this 
paper. Now using the general formula (1141) . we can com¬ 
pute the corresponding transition probabilities 

Po-y 2 = |Co—> 2 | 2 (15) 


taking various template of gravitational wave-forms that 
are likely to be generated in runaway Astronomical 
events. 

We start with the simple scenario of periodic GW 
generically written as 

hjk (t) = 2/ 0 cos fit (e x rf k + e+a-j k ) (16) 

where the amplitude varies sinusoidally with a single fre¬ 
quency U. In this limiting case of an exactly monochro¬ 
matic wave the temporal duration of the signal is infinite 
and we get for the transition probability 

Po ->2 = (tt/oU ) 2 (e x 2 + A 2 e+ 2 ) 

x [<5 (2 h7 + fl) — S (2m — fl)] 2 (17) 

where 

(18) 

is a dimensionless parameter carring the NC signature. 
Here we have used the stringent upper-bound^ |0| ~ 
4 x 10~ 4O m 2 for spatial noncommutativity and for ref¬ 
erence mass and frequency used values appropriate for 
fundamental phonon modes of a bar detector— which are 


formally identical to the NC harmonic oscillator system 
considered here. 

Consider periodic GW signal coming from a binary sys¬ 
tem (with quasi-circular orbit) being received by some 
earth-bound detector; if the orbit of the binary sys¬ 
tem is edge-on with respect to us, then we receive the 
-l-polarization of the radiation only£ , i.e., (e Xl £+) = 
(0,1), and in this case (fTTl) shows that the transition 
probability will scale quadratically with the dimension¬ 
less parameter A characterized by spatial noncommuta¬ 
tivity. Therefore such a transition will be driven by the 
combined perturbative effect of GW as well as spatial 
noncommutativity^ and will occur only if the space has 
a NC structure. In other words, a quantum mechanical 
transition induced by the linearly polarized GW from a 
binary system with its orbital plane lying parallel to our 
line of sight can be an effective test of the noncommuta- 
tive structure of space. Note that the angular frequency 
H of the quadrupole radiation is twice the angular fre¬ 
quency of rotation of the source^, thus (fTTl) also tells us 
that for transition to occur we need to have a harmonic 
oscillator with natural frequency m that matches with that 
of the source. Highly accurate X-ray/radio-astronomical 
measurement of the frequency of orbital rotation of bi¬ 
nary Pulsars can be used to pin-point the natural fre¬ 
quency of the harmonic oscillator required here. Since 
response of most operating bar detectors^ at present 
are centred around 0.9kHz, we can alternatively assess 
the target source-type by estimating the orbital radius of 
the source binaries which can emit GW of a suitable fre¬ 
quency to enable the said transition. Using Kepler’s law 
to express the orbital radius in terms of the angular fre¬ 
quency of rotation (= natural frequency of the harmonic 

oscillator in the present case) and total mass of the bi- 

1 

nary system, we find R = j - - '' tL-dlzl | 3 ss 23km where 

TOi = m 2 = 1.4 Mq, which is the lower mass-limit for 
neutron stars in terms of the Solar mass Mg, have been 
used as reference. Such small orbital separation can only 
be reached by highly compact binaries consisting of neu¬ 
tron stars and black holes^. 

For completeness let us mention that the amplitude / 0 
in eq. (flfil) for quadrupole radiation^ of GW from an in¬ 
spiral binary with masses mi, ni 2 and quasi-circular orbit 
is 



(19) 


where M c = ( mim2 l ? i j s the c hi r p mass of the binary 

(mi+m 2 ) 5 

and r is the distance of the binary from the Earth. 

From another binary system similar to the one consid¬ 
ered above, but with its orbital plane perpendicular to 
our line of sight, both the + and x polarization of the 
radiation will reach the detector with equal amplitude 
and consequently we will have a source for circularly po¬ 
larized GW signal that can be generically written as 

hjk ( t) = 2/ 0 [e x ( t) a) k + e+ ( t) cr%\ (20) 
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with e_|_ (t) = cos fit and £ x (t) = sin fit and amplitude 
given by eq. m- The transition probability in this case 
is 


P 0 ^2 = (tt/o^) 2 {(l+A)5(2w + fi)} 2 


+ {(1- A) <5 (2w - fl)} 2 


( 21 ) 


First let us note that as in the previous case, the tran¬ 
sition occurs when natural frequency of the harmonic 
oscillator resonates with that of the source of the GW. 
Also notable is that here the transition probability has 
terms both linear and quadratic in the dimensionless NC 
parameter A. However eq. m shows that for phonon 
modes in a bar detector which are the realization of the 
NC HO system considered in this paper, A is of the or¬ 
der of unity£2, so we cannot drop the quadratic term 
even though we started with a theory to first order in 
the NC parameter. Significantly, eq. m shows a non¬ 
zero transition probability for A = 0, i.e. if our space 
has commutative structure. Thus a transition induced 
by circularly polarized GW from a binary system cannot 
be used as a deterministic probe for spatial noncommu¬ 
tativity. This feature lies with the earlier case of linearly 
polarized GW signals only. 

In the last stable orbit of an inspiraling neutron star 
or black hole binary or during its merging and final ring- 
down, the system can liberate large amount of energy in 
GWs within a very short duration 10 _3 sec < r g < lsec. 
Such signals are referred to as GW bursts. Supernova ex¬ 
plosions and stellar gravitational collapse are other candi¬ 
date generators. Since bursts originate from violent and 
explosive astrophysical phenomena, their waveform can¬ 
not be accurately predicted and only be crudely modelled 
as 


hjk (t) = 2f 0 g it) (e x a) k + £+cr 3 fc ) (22) 

where to be generic we have kept both components of the 
linear polarization. Here g (t) is a smooth function which 
goes to zero rather fast for |t| > r g . A convenient choice 
is a function peaked at t = 0 with g (0) = 0( 1) so that 
| hjk (f) | ~ O (/o) near the peak. So we take a simple 
Gaussian 


where g (fl) = \fKT g e V 2 ) is the amplitude of the 
Fourier mode at frequency fl. 

Using ea. (124l) in the general formula for transition am¬ 
plitude (fT4l) and further employing ea. dlhll . we find the 
probability for transition from the ground state to the 
second excited state induced by a GW burst is 

P 0 ^2 = (2 v ^/ 0 n7T g ) 2 e- 2ro2 U 2 ( £x 2 + A 2 £+ 2 ) (25) 

where the GW Fourier mode with twice the natural fre¬ 
quency of the harmonic oscillator (the detector in our 
consideration) gets picked up. Since the burst signal du¬ 
ration T g ~ 10 -2 — 10 _3 sec, the maximum frequency in 
the Fourier spectrum can be H max /27r ~ 0.1 —1kHz which 
partially overlaps with the sensitive bandpass for the bar- 
detectors. Actually the lower limit is in the bandpass of 
LIGO which, though, does not represent our harmonic 
oscillator system as closely as a bar-detector, still works 
in the long wavelength-low velocity limit considered here. 

From ea. (l?5l) we again see that the + polarization of 
the GW burst can only induce a transition if the space 
has a NC structure. As mentioned earlier the polariza¬ 
tion state of a GW signal from a given source can be 
anticipated since it depends largely on the orientation of 
the source which can be determined by observing its elec¬ 
tromagnetic radiation and the detector geometry which 
is at our disposal. So detecting a QM transition induced 
by a GW burst from an appropriate source can serve as 
a probe of the spatial noncommutativity. 

A slightly more realistic waveform can be realized if 
the Gaussian function in ea. (l23D is modulated by some 
frequency Ho/27r, resulting in a sine-Gaussian amplitude 

g(t) = e -t /r s sinS2 0 t . (26) 

Then the signal in eci. fliHl) admits the Fourier modes 


9(0) = 


iy/nfoT g 


e -(n-o 0 ) 2 r 3 2 /4 _ e -(n+n 0 ) 2 r 2 /4 


centered around fl = flo in the frequency space, 
corresponding transition probability is 


Po-s-2 — e 


_ (B -(2ro-n 0 ) 2 r 2 /4 _ e -(2ro+n 0 ) 2 T 2 /4 


(27) 

The 



(23) 


x (v^/oWTg) 2 ( £x 2 + A 2 £ + 2 ) (28) 


Owing to its small temporal duration the burst have 
a continuum spectrum of frequency over a broad range 
upto / m ax ~ 1/Vg whereas the detector is sensitive only 
to a certain frequency window and blind beyond it. If 
the sensitive band-width is small compared to the typi¬ 
cal variation scale of the signal in the frequency space, the 
crude choice in ea. (!23D instead of a precise waveform is 
good enough. In terms of the Fourier decomposed modes 
the GW burst can thus be modelled as 



where, as usual, only the appropriate Fourier mode is 
picked up by the harmonic oscillator. Two limiting cases 
can be immediately identified in ea. (E51) . If the operat¬ 
ing frequency of the detector is low, say, in the sub-Hz 
bandpass, the two exponential terms are of same size and 
nearly cancel each other, reducing the transition proba¬ 
bility. The other extreme is when 2m — flo = A m with 
Ap << 1. Then the first exponential term is sizeable 
while the second term is negligible in comparison and we 
get a transition amplitude 

Po—s-2 ~ 


g{fl)e~ int dfl (24) 


e-(A^r>/2 foWTg f ( £x 2 + A 2 £ + 2 ). (29) 
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Note that since r g ~ 27r/fi max , the quantity (Aror g ) 2 ~ 
(27rA w/Q max ) . For bar detectors Aw is at least within 
the sensitive bandpass^, so it cannot be more than a few 
Hz, whereas fl max is in the kHz range for GW burst of 
a few millisecond duration. This ensures that the Gaus¬ 
sian factor in ea. (l29l) suppresses the transition probabil¬ 
ity only marginally. Interestingly, the Gaussian factor 
in ea. (1251) is much stronger compared to that in ea. (l29D . 
since ror 9 >> A bt 9 . Thus for the more realistic sine- 
Gaussian template, the transition probability is actually 
higher. 

In conclusion we would like to convey that the con¬ 
siderations in the present paper suggest that the joint 
operation of various resonant detector groups like ALLE¬ 


GRO, AURIGA, EXPLORER, NAUTILUS and NIOBE 
around the world in IGEC (InternationalGravitational 
Event Collaboration)^ - — may possess the potential to 
establish the possible existance of a granular structure of 
our space as a by-product in the event of a direct detec¬ 
tion of GW and therefore must be continued. 
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